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We investigate the fidelity of the quantum state transfer (QST) of two qubits by means of an arbi¬ 
trary spin-1 network, on a lattice of any dimensionality. Under the assumptions that the network 
Hamiltonian preserves the magnetization and that a fully polarized initial state is taken for the lat¬ 
tice, we obtain a general formula for the average fidelity of the two qubits QST, linking it to the one- 
and two-particle transfer amplitudes of the spin-excitations among the sites of the lattice. We then 
apply this formalism to a ID spin chain with XX-Heisenberg type nearest-neighbour interactions 
adopting a protocol that is a generalization of the single qubit one proposed in Ref. [Phys. Rev. 

A 87, 062309 (2013)]. We find that a high-quality two qubit QST can be achieved provided one 
can control the local fields at sites near the sender and receiver. Under such conditions, we obtain 
an almost perfect transfer in a time that scales either linearly or, depending on the spin number, 
quadratically with the length of the chain. 


I. INTRODUCTION 

The capability of faithfully transferring information 
from one location to another is one of the main driv¬ 
ing factors of the modern technological progress. As far 
as classical information is concerned, there is no limit, 
at least in principle, to reproduce an exact copy of the 
original message; therefore the information transfer has 
to face somehow minor problems than those faced in the 
quantum realm. There, the no-cloning theorem [1] ex¬ 
plicitly prohibits to make an exact copy of the quantum 
state on which the quantum information has been coded 
in. This has stimulated, over the last past decades, a 
large body of works on how to efficiently achieve Quan¬ 
tum State Transfer (QST). 

For short-haul transfers of the quantum state of a sin¬ 
gle qubit (1-QST), the use of spin-| chains, initially 
proposed in Ref. [2], has been largely investigated (see 
Refs. O [4] and references therein, and Ref. [5] for an 
implementation with a cavity array). Protocols based 
on time-dependent couplings mil], fully engineered in¬ 
teractions [HE], ballistic transfer [nHni, Rabi-like os¬ 
cillations USES], just to name a few, have been shown 
to achieve high fidelity 1-QST, in addition to some addi¬ 
tional tasks like routing of the quantum information to 
an on-demand location on a spin graph [29ti3T] . 

Recently, the same effort is being devoted to the case 
of multiqubit QST (n—QST), in which the state aimed 
at being transferred is made of n > 1 qubits. In many 
cases, the adopted strategies consist of extensions of 1- 
QST protocols and, as a consequence, the drawbacks and 
inconveniences they already presented for the 1-QST are, 
to some extent, even more amplified when it comes to the 
n-QST case. For example, the multi-rail scheme [32l [33] 
requires the use of several quantum spin-^ chains and a 


complex encoding and decoding scheme of the quantum 
states; employing linear chains made of spins of higher 
dimensionality reduces the number of chains to one, but 
still requires a repeated measurement process with con¬ 
secutive single site operations [34]; the fully-engineered 
chain (eventually combined with the ballistic or Rabi- 
like mechanism), as well as the uniformly coupled chain 
with specific conditions on its length, needs conditional 
quantum gates to be performed on the recipients of the 
quantum state [35f[37| . Therefore, simpler many qubits 
QST schemes would be quite appealing. 

In the present paper, we adopt a minimal engineer¬ 
ing and intervention point of view, looking for a 2- 
QST protocol that does not need demanding operations 
to be performed, neither in the form of external end- 
operations on the spins nor to engineer the spin cou¬ 
plings. Experimentally-friendly 2-QST schemes are in¬ 
teresting in view of the fact that both the full modulation 
of the couplings may be unattainable (depending on the 
physical system meant to perform the QST) and that 
quantum operations, such as measurements and gates, 
are prone to errors which, in a realistic set-up, may fatally 
degrade the efficiency of the protocol. In addition, as the 
exchange of quantum information is meant to occur, for 
instance, between quantum processors, it is quite natu¬ 
ral that QST of more than a single qubit has to be faced 
in order to fully exploit the potentialities of quantum 
computation. Needless to say, also other fields relying 
on quantum information processing, such as cryptogra¬ 
phy and dense coding would widely benefit from efficient 
n-QST protocols [38] . 

The paper is organized as follows: in Sec.[ll| we obtain 
a general expression for the average fidelity of the quan¬ 
tum state transfer of two qubits coupled to an arbitrary 
total angular momentum conserving graph of spin-| ini- 
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tialized in the fully polarized state; in Sec. |m]a specific 
one-dimensional instance of such a graph is presented and 
it is shown that, by means of strong local magnetic fields 
on the so-called harrier qubits [23], high-quality 2-QST 
can be achieved. Finally, in Sec. [TV| some concluding re¬ 
marks are reported together with a discussion on possible 
extensions of our idea. 

II. FIDELITY FOR A CLASS OF SPIN-| 
HAMILTONIANS 

In this Section we derive a general expression for the 
average fidelity of a 2-qubit quantum state transfer from 
a pair of senders to a pair of receivers, residing respec¬ 
tively on sites S = {si, 52 } and IZ = {ri, r 2 } of a lattice K 
of arbitrary dimensionality. The only constraints we as¬ 
sume to be satisfied by the spin dynamics on 1C are 1.) the 
conservation of the total magnetization = YlneK: 
along some axes a (which we assume hereafter to be the 
quantization axes z) and 2.) the initialization of all the 
spins of 1C but S into a fully polarized state along 2 ;. 

The most general Hamiltonian, allowing up to two- 
body interactions, for spin-^ particles is given by 

( 1 ) 

ij£lC a /3 

where jS = {0,x,^,^} with 5'^ = |. Because of 
the conservation rule implied by = 0, Eq. 0 

can be decomposed into a direct sum over all subspaces 
with fixed 2 ;-component of the angular momentum, H = 
0^^ Sz = - Y: •••5 T- Without loss of gen¬ 

erality we re-scale the labelling of the angular momentum 
sectors by the number n of spins flipped in each sector, 
that is Sz = —Y with n = 0,1, 2,..., N. The Hilbert 

space dimension of the n-th sector is clearly 2^^). 


Our goal is to transfer the quantum state 
of two qubits located at sites S and given by 
|' 0 (O))^ =(a | 00 )+/3 | 01)+7 | 10 )- 1 -^ | 11 ) to the receivers 
spin, located at sites IZ. The rest of the lattice, embod¬ 
ied by the quantum channel F and the receivers IZ^ is 
initialized in the state \TTZ) = |0)^-. The evolution 

of the overall state |^(t)) is given by 

= 100)5 |^ 7 ^) + 

+e -*^^*(/3 101)5+7110)5) |r 7 e) + 

111)5 |^7^), ( 2 ) 

where the Hamiltonian has been restricted to the sub¬ 
spaces n = 0 , 1 , 2 , respectively, by taking into account 
the invariant sector of the Hilbert space to which each 
component of the state vector pertains. 

By tracing out all of the spins but the re¬ 
ceivers, one obtains the state of the latter, P 7 ^(t) = 
TrK:- 7 ^ (|^(t))(^(t)|). The fidelity between the state 
transferred to the receivers and the state encoded ini¬ 
tially on the senders is given by [39] 

Fmmmis >= simipnit) i^(o))5. (s) 

The quality of a QST protocol, however, cannot be 
simply evaluated by considering the fidelity of the trans¬ 
fer of a single, specific input state; in fact, a more appro¬ 
priate figure of merit is given by the average QST-fidelity 
F{t) obtained by averaging over all possible input states. 

After a lengthy but straightforward calculation, full 
details are reported in [40] we obtain the average fidelity 
F{t) for the 2 -QST with the constraints of a lattice 1C 
described by a total 2 ;-magnetization conserving Hamil¬ 
tonian and provided the fully polarized initial state is 
taken for F and IZ : 


F{t) = 


4 54 


fri _j_ Fr2 _ fr2 ( i’ri 

Js-\_'Js2 ' S2\J S\ ) 


(i/i'/P+i/r/P) (i/i'a P+i/r/P) [dsiz] + 


/ 71,171^71 \ ^ 1 

-5J 1- E +5jb» [(/;;+o wg)*]E [(/.“.)•«.+(/;)•«.]. 

\ n<m=l / n=l 


(4) 


where /^= {m\ \n) and g^m— (^'^1 \nm) are 

the single- and two-particle transfer amplitudes from sites 
n ^ m and {nm} {rs}, respectively. 

Eq. @ plays the same role for the 2 -QST of the cele¬ 
brated average fidelity expression given in Ref. [ 2 ] for the 
1 -QST. 

Notwithstanding the lengthy expression for the aver¬ 
age fidelity, in the presence of further symmetries and 


specific Hamiltonians intended to implement the 2 -QST 
protocol, Eq. 0 can be considerably simplified. In the 
next Section we give an instance of such a procedure and, 
at the same time, we propose a model that accomplishes 
a high-quality 2 -QST. 
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III. THE MODEL AND THE PROTOCOL 

The results for the 2-QST scheme we propose in this 
Section have to be compared with the average fidelity 
(hereafter called fidelity) we would obtain by means of 
local operations and classical communication (LOCC) 
or by means of universal quantum cloning machines 
(UQCM). The use of these channels yields what is con¬ 
ventionally dubbed as classical fidelity and amounts to, 
respectively, ^LOCC = ^|q; HI] and ihjQCM = 2(1m) @2’ 
where d is the Hilbert-space dimension of the state aimed 
to be transferred. For the case of 2 qubits we have 
d=A and, therefore, our protocol outperforms the clas¬ 
sical ones if we obtain a fidelity higher than | (or | if 
optimal cloning is not available for the system at hand). 

The lattice 1C we will consider is a ID spin-^ open chain 
and the Hamiltonian is taken of the AX-Heisenberg type 
with nearest-neighbor interactions only, and a magnetic 
field along the z-axis on the 3rd and {N — 2)th spin, 
playing the role of the ‘barrier’ qubits, separating the S 
and IZ pairs from the rest of the channel: 


N-l 

if = - y] J((crfcrf +1 + o-fo-f+i) +h{al+ a%_ 2 ) , (5) 

1 = 1 

where = 25'^ {a = x, y, z) are the usual Pauli matri¬ 
ces. 

We aim to achieve the transfer of an arbitrary 2-qubit 
state residing on the sender spins 5, located at sites 1 and 
2, |'0(O))^2 |00) +/^ +7 111 ) 5 1® receiver 

spins 7^, residing at the other end of the chain, IZ = 
X — 1, X, as depicted in Fig. 


h 
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FIG. 1: (Color online) The spin graph JC by means of which 
we aim at achieving high-quality 2-QST via a Rabi-like mech¬ 
anism between the two ends of the spin chain. The quantum 
state to be sent is encoded on the spins S = {1?2} (blue), 
whereas the receivers are located at X = {X — 1,X} (red). 
Barrier qubits, residing at sites n — {3,X — 2}, on which a 
strong magnetic field h is applied are colored in green. 


Eq. (§ can be mapped to a tight-binding spinless 
fermion model via the Jordan-Wigner transformation [43] 

AT-l 

-H" = - E 4ci+l+Cic|+i+/l (4 c3+c]v_2CJV-2) , (6) 

where we have taken as our energy and inverse time unit 
the exchange energy J, that we consider to be site inde¬ 
pendent. 

Because of the quadratic nature of the Hamiltonian, 
the single particle spectrum is sufficient to describe the 
full dynamics. Denoting by Sk and \ek) the k-th en¬ 
ergy eigenvalue and its corresponding eigenvector, the 
full Hamiltonian operator acting on a 2^ dimensional 
Hilbert space, is easily decomposed into a direct sum 
over all particle number-conserving invariant subspaces 
H = 04=1 Hn, where 


N 

Hn= E l{ 0 }) ({0}|Cfe„-Cfe2Cfci, ( 7 ) 

k\<k2<.-.<kn = l 


with |{0}) being the fermion vacuum. Each therefore 
can be constructed quite simply once the single-particle 
spectrum is known. Notice that the specific ordering of 
the kCs in the sum of Eq. 0 is taken in such a way 
that unwanted phase factors do not arise when mapping 
back into spin operators via the inverse Jordan-Wigner 
transformation. 

Therefore, in order to evaluate |4^(t)) as given 
by Eq. 0 . we need the spectral resolution of 
iJi, given by the eigenvalues and eigenvectors of 
the following X x X tri-diagonal matrix: Tij = 

which is easily diagonal- 
izable, at least numerically. Notice that a uniform mag¬ 
netic field along the ^-direction has no influence on what 
follows as it corresponds to adding a term proportional 
to the Identity in Eq. Hence the eigenvectors do not 


change, whereas the uniform shift experienced by all of 
the eigenvalues is cancelled out in the time evolution of 
the fidelity, which, as we will show below, only depends 
on energy differences. 

Key to our aim is the presence of eigenstates that 
are at the same time strongly localized on both the 
sender and the receiver spins. That is, by expand¬ 
ing the Hamiltonian eigenvectors in the position basis 

kfe) = Z:4=i 1^) (where \n) = |0i....0„_il„0„+i...0jv) 

describes a state with a single spin flipped at position 
n), a prerequisite for Rabi-like oscillations based 2-QST 
protocol to correctly work is that there exists a certain 
(small) number of eigenstates \ek) for which a^n is non- 
negligible only for n = 1,2,X—1,X. We find this re¬ 
quirement of edge quadri-loealization to be fulfilled for 
spin chain of lengths X7^3n—1, where n > 3 G IN, when 
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strong magnetic fields are applied on the barrier qubits. 
The condition n > 3 is due to the fact that the minimum 
length of a spin chain allowing for 2 senders, 2 receivers, 
and 2 barrier qubits is N > 6 . 

By writing Hi = k/c)(^/c| with the eigenvalues 

taken in increasing order, the localized states are labelled 
hyk = {Q[7V,3]-1,Q[7V,3],7V-Q[7V,3]-1,7V-Q[7V,3]}, 
where Q[a,b] denotes the quotient. In the following we 
will refer to these states by {^i}, i = 1,2, 3,4. In Fig. 
an instance of such localized structure of the eigenstates 
is given for N = 46 and h = lOOJ. We observe that 
there are 4 eigenstates, labelled by k = 14,15,30,31, 
that are quadri-localized on the edges, i.e. at sites n = 
1,2,45,46. Besides these four eigenstates, another two 
are bi-localized on the barrier qubits at sites n = 3,44; 
whereas the remaining ones are extended states with neg¬ 
ligible amplitudes on the senders, the receivers, and the 
barriers. As a consequence, the contribution to the dy¬ 
namics of |T(0)) of these extended states is negligible up 
to order 0{h~^). In the case of A" = 3n — 1, on the 
other hand, two additional extended eigenstates appear, 
labelled by k = {Q[A, 3] +1, A —(5[A, 3] —2} with a non- 
negligible value of a^n for n = 1, 2, A — 1, A, as shown in 
Fig. for the case of A = 50. We will refer to the latter as 
extended edge-localized states. As a consequence, there 
are more eigenstates taking part in the time evolution of 
the initial state and, although, high-quality 2 -QST is still 
attainable, the clear-cut analysis we will give below is, to 
some extent, complicated by the their presence. There¬ 
fore, in the following, we will first consider spin chains of 
length A 7 ^ 3n — 1 


channel F (due to the energy mismatch with the con¬ 
necting spins at sites 3 and A — 3), their energies and 
eigenstates can be approximated by Ist-order degenerate 
perturbation theory and read 
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Qi 



— (^qin\'kk) , |ag^^|— ^5 

n={l,2,Ar-l,Ar} 


( 8 ) 


where i = 1 , 2 ,3,4. Notice that the coeffi¬ 

cients of the eigenvectors obey the parity relation 
Cikn={—^Y^^CikN+i-n becausc of the mirror symmetry 
of the model [44] . 


Exploiting the quadratic nature of Eq. (|^, we can 
reduce the two-particles transfer amplitude to single¬ 
particle ones by means of the relation given in Ref. [36l 

mm 


Pr y’r+l , _ f 

J n J n ' ' ' J1 


s — 1 


9 


rs 

nm 


-pr 

/n+1 


f^+i 


J m ' ' ' J m 


( 9 ) 


A. Rabi-like 2-QST 

Because the quadri-localized states come as a result of 
the small effective coupling of S and IZ to the quantum 


Moreover, mirror symmetry millT] implies |/[| = 
and perturbation theory allows to retain only the tran¬ 
sition amplitudes between the senders and the receivers. 
Working out all these simplifications the average fidelity 
given by Eq. Q reduces to the approximate expression 


= i + [(/"-■)"] + ^l/i 
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nN-l* nN nN-1 

J 2 


N-l\2 I ^ I ^Ar|2 

+ ^l/i I 


( 10 ) 


r 


which we remind to be correct up to order O [h ^). 


because of Yln=i \fr\‘^ — for alH G 1C coming from 


Eq. ( [Tq| ) will be the starting point of the following anal- the conservation of M^. Although /i (t), /i (t), and 

are complex-valued functions of time, Fa{t) is 
a real-valued bounded function, which, taking Re[/J(t)] 
and Im[/J(t)] as independent, becomes a function of six 
real-valued bounded functions. Therefore, standard La- 
grangian multiplier methods can be applied in order to 
search for the absolute maximum of Fa{t) within the 
boundaries given by Eqs. (11). It turns out that the max¬ 
imum of Fa{t) is given by the conditions Re[f^~^{t)] = 1 


ysis in which we will evaluate both the maximum achiev¬ 
able fidelity and the optimal transfer time. To start with, 
let us notice that Fa{t) only depends on the three com¬ 
plex variables and which obey the con¬ 

straints 




( 11 ) 
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FIG. 2: (Color online) Density plot of \akn\ versus k and n 
for N = 46 (left) and iV = 50 (right) with h = lOOJ. In the 
case of N ^ 3n — 1 (left panel) it is clearly shown that the 
spatial distribution of the eigenvectors exhibits the presence 
of 4 quadri-localized energy eigenstates; whereas, in the case 
of N = 3n — 1 (right panel), two more eigenstates appear 
with non negligible amplitudes also on sites S and IZ. Their 
presence gives rise to a more complicated dynamical behavior 
of the hdelity because of the larger number of degrees of free¬ 
dom (and frequencies) effectively entering the time evolution 
of |T(0)). 


and Re[fl^{t)] = (t)] = /m[/J(t)] = 0 (the lat¬ 

ter following from the former due to the conservation 
of M^) and amounts to Fa{t) = || :^ 0.97. We found 
that Fa{t) does not achieve the maximum possible value 
of 1 just because it is an approximate expression for the 
fidelity: in fact, if the values obtained above for the tran¬ 
sition amplitudes are used in the exact expression of the 
average fidelity given by Eq. Q, we obtain F(t) = 1. 

The next step is to find the time t* at which the transi¬ 
tion amplitude reaches the optimal values for the 2-QST. 
To do this, we can maximize the function Re[fi~^{t)] (or 


equivalently, due to mirror-symmetry, i?e[/^(t)]) whose 
time evolution is generated by the Hamiltonian given in 
Eq. (§. This will fix the transfer time t*. 

Although Fa{t) is an highly oscillating function be¬ 
cause of the presence of many frequencies in the transi¬ 
tion amplitudes, it is possible to find the transfer time t* 
of the 2-QST protocol in a relatively simple way, as we 
will outline in detail in the following, for the case of even 
N^3n—1. 

By exploiting Eqs. <© and by means of elementary 
trigonometric identities, the term Re[fi~^] can be ex¬ 
pressed as 

N 

fc=i 

4 

~ ( 12 ) 

i=l 

which, for even N, becomes 

- (_l)Mod[iV,3]+l 

X (sin cJq t cos cJq t sin t cos t + 
-hcoscjQ t sin cJq t cos sin , (13) 

where 


CJn = 


CJ, = 


UJ 


14 ■ 


■ UJ. 


23 


u;i4 - uj. 


23 






(14) 


and Mod [a, b] is the modulus function. 

Having defined in such a way the frequencies that 
enter the dynamics of Re[f^~^]^ it turns out that 
, which immediately sets a time scale 
for the 2-QST. In fact, as we can focus only 

on the first summand of the RHS in Eq. namely 
sinujQ t cos sin ujit. On the same footing, the next 
time scale is given by uj^ , which implies that the maxi¬ 
mum of Re[f^~^] (t) has to be found in the neighborhood 
of ti = -M, which, hence, approximately gives the opti- 

mal time t*. As a result, the transfer time t* is found by 
solving 


t* 


max [sin uJq t cos t] 
min [sinco’^tcosco’^t] 


if Mod[N,3] is odd 
if Mod[N,3] is even 


(15) 

and by choosing the solution closest to the time ti. A 
graphical representation of these timescales is given in 
Fig-i which also allows us to put forward a simple phys¬ 
ical interpretation. In fact, from the left panel of Eig. 
sin uj it can be seen as an approximate envelope for the 
fidelity of the transfer process, whereas uJq 2 J gives 
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the time-scale for the (very rapid) bouncing of the exci¬ 
tation back and forth between the two receivers. These 
oscillations occur because of the direct coupling between 
the two receiver spins and could be eliminated if this 
coupling is switched off. 

This bouncing occurs many times (see the right panel 
of Fig. 1^ before the excitation slowly leaves these two 
sites to go back towards the senders. In Fig. we 
compare the approximate value of the fidelity given by 
Eq. ( p^ with the exact result of Eq. 0 : it is shown 
that F\t)>Fa{t) attaining values as high as 0.999. We 
checked, up to computational accessibility, that this holds 
true for every N. 


proximate value of the optimal time for the excita¬ 
tion transfer from S to IZ^ can be obtained by using 
standard degenerate time-independent perturbation the¬ 
ory to evaluate the relevant energy eigenvalues. Here 
we report the energy corrections for the dynamically 
relevant states given in Eq. (§ up to the first order 
in h~^. Denoting by with xi>X 2 >x^^ the solu¬ 
tions of —and defining the param¬ 
eters = ^+^ 1 , 2 , OL = — 1 , and 7 = 

(2 (<aq 2 +/^i, 2 +l)) ^, we obtain 

^qi ~ 5 ^q2 ~ 5 ^^3 ~ ^^2 5 ^^4 ~ ^^2 ’ 




FIG. 3: (Color online) (left) Approximated average fidelity 
Fa(t), and, for comparison, |sinct;7t| shown as a function of 
time. Notice that Fa is so rapidly oscillating that it appears 
to fill the entire red region in the plot. The sinusoidal function 
approximately gives an envelop of the fidelity with the same 
periodicity, which easily allows us to identify several reading 
windows for the 2-QST. (right) Full expression for the av¬ 
erage fidelity Fa{t) (red), sintcost (blue) and sincj^t 
(green) shown as a function of time around the optimal trans¬ 
fer time ti, denoted by the vertical dashed line. Notice that 
|sinct;7t| — 1 and that the fidelity attains very high values 
near ti which recur several times, with a frequency of order 
of allowing again for several reading windows. 

The Rabi-like half-oscillation time ti, giving an ap- 


where 


a7 = 2 


Zl,2 ■ 




N-6 

A^-5 E 

k=l 


T'1.2 


((1 ± COS kn) sin 


Zl,2 


2coSj^ 


/ 


(17) 

we obtain the ap- 
_A++A 2 --A+)]. We 


Finally, using Eqs. (14) and ( |16[ ) 
proximate transfer time [2 
also find that ti scales quadratically with the mag¬ 
netic field intensity and that it is independent of TV, 
ti = |h^ + (—as reported in Fig. 
for the case of N=30. Notice that in 1 -QST schemes in 
which the magnetic field is applied directly on the sender 
and the receiver, the transfer time scales exponentially 
both with the length of the chain N and with the mag¬ 
netic field’s intensity h [191 [ 22 ]. 

A similar procedure for odd N^3n — 1 yields for the 
transfer amplitude of Eq. 

X (^cosuJq t sin uj^t sin cos uj ft 

+ sin uo^t cos t cos sin uo ft) ,(18) 

and the optimal transfer time t* can be found via a 
double step, with a recipe similar to the one discussed 
above. Eirst we determine ^ 2 , defined as the solution 
of sinco’^t=(—which is closest to = and 

then is given by the solution of cosci;^t=l which is 
closest to ^ 2 - Notice that this is not different from the 
previous procedure; indeed, there are many ways to re¬ 
cast Eq. (12) by combining the energies 5 ^., and the fact 
that we obtained an apparently different method for t* in 
the Eqs. (13) and ([l^, for even and odd AT, respectively. 


is due to the fact that we kept unchanged the definitions 
of the frequencies given by Eq. (14) to avoid a confusing 
re-labelling. 


B. Quasi Rabi-like 2-QST 

Let us now go back to the case of spin chains of length 
AT=3n—1, that was left out of the previous analysis. In 
this case, two additional edge-localized extended states 
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FIG. 4: (Color online) (left) Exact average fidelity F{t) (yel¬ 
low) and approximate fidelity Fa{t) (red) versus t around 
ti. Notice that F(t) > Fa{t) is always fulfilled and that, 
at the optimal time almost perfect quantum state transfer 
is achieved, with F{t) ~ 0.999. (right) Transfer time t* ver¬ 
sus h for iV = 30: the function f{h) = f (dashed line) fits 
perfectly with the numerical data for t *. The range of values 
of the magnetic field h is such that Fa{t) > 0.95. 


are found, whose presence hinders the clear Rabi-like os¬ 
cillations of the excitation between S and 7^, exhibited in 
the case of 1 by the sinusoidal function sinco’ft 

and reported in the previous subsection and in Fig. 
Nevertheless, since the number of eigenstates (eigenener- 


gies) to be included in the sum given in Eq. 12 increases 
just by two, an analysis similar to the one performed 
above can still be carried out. We dub the transmission 
process in this case as quasi Rabi-like 2-QST. 

In fact, the approximated expression for Re[f^~^{t)] 
given in Eq. (12) now becomes 


6 

i—\ 

= ±2 (ci cos cos 
— C 2 cos UJ 2 ^t cos + 

+ C 3 coscj^t coscj^t) , (19) 


with the ± sign holding for even (odd) N, respectively. 


and where we exploited the mirror-symmetry and used 
Ist-order degenerate perturbation theory relations 

1 3 

Cl — aiiaiTv-i — a4ia4Ar-i — ^ _ i ’ 

C2 = Ci2ia2N-l — (21) 

3 

Using the approximations —2 and 

- 0> yields 

Re[f^~^] ~ ±cos2t ( 4 c 3 sin — sin^ . (23) 

Once again, it is possible to identify two different pro¬ 
cesses: the slow quasi-Rabi-like oscillations of the exci¬ 
tation between S and 71, having a time scale ruled by 
and the fast oscillations of the excitation within IZ 
(S) triggered by J and described by the term cos 2 t. Al¬ 
though the additional states complicate somehow the ex¬ 
pression of Re[f^~^], they also provide a clear advan¬ 
tage as far as the transfer time is concerned. Indeed the 
time ti is now linear in the magnetic field and hence 
2 -QST occurs faster w.r.t. chains of length Ny>^3n—1. 
Let us remind that a similar phenomenon takes place in 
1 -QST protocols where the sender and the receiver are 
weakly coupled to the chain either because of smaller 
bond strengths [ 20 ] or because of a strong magnetic field 
acting on barrier qubits [23]. Indeed, in those cases the 
QST time for even- and odd-length chains too scales, 
respectively, quadratically and linearly with the pertur¬ 
bation’s intensity. 

In Fig.[^ we summarize the main result of the two pre¬ 
vious subsections, namely the possibility to transfer with 
high-fidelity an arbitrary quantum state of two qubits by 
means of a linear spin- ^ chain with strong magnetic fields 
on the barrier qubits. It is shown that for A^ 7 ^ 3 n—1, a fi¬ 
delity close to unity can be achieved regardless of N (pro¬ 
vided strong enough magnetic fields are applied at the 
barrier sites) although in a time that increases quadrat¬ 
ically with h. For A^=3n—1, instead, the quality of the 
transfer depends on whether N is divisible by 4 (higher 
Fa{F)) or not (lower Fa{t*)), but the transfer is achieved 
in a time that scales only linearly with both h and N. 
Nevertheless, the differences amongst all these cases fade 
away for where all curves collapse. 

Since the average fidelity is not identically one, one 
could imagine that there exist specific input states that 
are transferred with a relatively poor quality. This is not 
the case, and in order to dispel such a doubt, we evalu¬ 
ated also the worst case fidelity m and found that the 
minimum state-dependent fidelity, evaluated by means of 
Eq. remains close to the average one up to the fourth 
digit, i.e. Fmin - 0.999. 

To conclude this Section we stress that the probabil¬ 
ity to find the excitations inside the quantum channel 

r, evaluated by ^ 2^=3 (\fiWf + l/:?(i)l^)) is radically 










jrxlO^ 
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FIG. 5: (Color online) (left) Transfer time for the case 
of spin chain’s of length N=3n—1 as a function of h. The 
time t* increases linearly with h, although with slightly dif¬ 
ferent coefficients depending on wether N is divisible by 4 or 
not. In the figure the two cases are exemplified by N=32 and 
A/'=38, respectively, (right) Transfer time t* for the case of 
spin chain’s of length iV=3n—1 with magnetic field /i=4000 
as a function of N. The time increases linearly with N, with 
clearly distinct coefficients depending on the divisibility by 4 
of iV. 


different depending on whether Ny^3n—1 or not. In the 
former case it is of the order of 0{h~^) because F acts 
as a mere physical connector |T9l [20l |26] entering the dy¬ 
namics only virtually. On the contrary, when N=3n—1^ 
two extended states with a non-negligible overlap on the 
senders and the receivers come into play, meaning that 
the excitations can be actually found inside the channel 
F. As a consequence, the effect of a dissipative coupling 
of an environment eventually acting only on F has a neg¬ 
ligible influence only for the Rabi-like QST whereas for 
the quasi-Rabi-like one the quality could be severely de¬ 
graded, especially for short chains. For long chains, on 
the other hand, since the overlap with the extended state 
localized also on the edges scales as C3 r\j O(Ar-i), 
degrading effect becomes negligible. On the other hand, 
the presence of disorder in the couplings or in a mag¬ 
netic field acting on the quantum channel F should have 
a negligible influence too on the efficiency of the 2-QST 
protocol we are proposing here, regardless of N, as long 



FIG. 6: (Golor online) Maximum of the approximate fidelity 
Fa{t*) obtained by choosing /i=4000. The red points are the 
Rabi-like 2-QST and the maximum does not depend on A, 
the blue point are for N—3n—l and even N performs better 
than odd N. For A^l, both of the curves converge to the 
N^3n—1 case. 


as the spatial distribution of the eigenvectors depicted in 
Fig. is not significantly affected j48H5Q] . 


IV. CONCLUSIONS 


In this paper, we derived an expression for the aver¬ 
age fidelity of the quantum state transfer of two qubits 
through a spin-| chain, providing that the z-total angu¬ 
lar momentum is conserved and all the spins are initially 
aligned. This general expression, relating the average 
fidelity explicitly to one- and two-particle transition am¬ 
plitudes, may result useful in investigating the two qubit 
QST properties of a wide range of physical models dis¬ 
playing the above-mentioned characteristics. 

Furthermore, we discussed a specific case, obtained by 
extending a Rabi-like protocol, widely adopted for QST 
of single qubits, to the non-trivial case of the QST of 
two qubits, where the senders and receivers pairs are lo¬ 
cated at each end of a one dimensional spin-1 chain with 
W-Heisenberg type nearest-neighbor interactions. The 
presence of strong magnetic fields on the two sites clos¬ 
est to the sender and receiver (barrier qubits), allowed 
us to obtain a faithfully transfer for an arbitrary two- 
qubits quantum state. We characterized the quality of 
the 2-QST by using first-order degenerate perturbation 
theory; thus providing, apart from a clear-cut physical 
interpretation of the multi-excitation dynamics yielding 
high-quality 2-QST, also an approximate analytical ex¬ 
pression for the time transfer. The latter is found to 
increase linearly or quadratically with the magnetic field 
intensity h depending on the spin chain length. More¬ 
over, we have verified that the worst-case fidelity of the 
QST remains almost unchanged with respect to the av¬ 
erage one, i.e., 0.999. 






9 


Two final comments are in order. Since it is straight¬ 
forward that the 4-dimensional Hilbert space of the 
senders can be also employed to encode an n < 4 qudit 
to be transferred to the receivers, the scheme we propose 
can be adapted to the transfer of qutrits (qubits) encoded 
in arbitrary three (two) orthogonal quantum states of 
the senders. Finally, the question if the barrier scheme 
discussed here is useful in order to perform QST of an 
arbitrary number of qubits will be left to future investi¬ 
gations. 
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